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SECTION 1 


INTRODUCTION 


1.1 BACKGROUND AND OBJECTIVES 


The detonation of one or more nuclear weapons on or near the 
Earth's surface generate a lofted cloud composed of dust and other 


debris. 
The lofted mass is a potential threat to 


e reentry vehicles which must fly through clouds from prior 
detonations (fratricide) 


@ ascent vehieles which must flyout soon after being attacked 


@ radar and communications which must operate in a 
dust/debris environment 


The objectives of this report are to provide interested per- 


sons with 
1. an overview of dust/debris cloud phenomenology 


2. examples of cloud characteristics 


1.2 BASIC CONCEPTS 


After a near-surface nuclear burst, the cloud characteristics 
of primary interest are the concentration* and particle size dis- 
tribution of the condensed phase mass in the cloud as a function 


of position and time. 
The condensed phase mass in the cloud can be composed of: 


e earth matertal which has been entrained and lofted by 
the ascending winds which develop. (Earth or ground 
material include soil/rock, water, vegetation, as well 
as structures near the burst point.) 


@ weapon material which has been vaporized and recondenses 
during the early phases of the nuclear cloud development. 


e water/ice particles which condense from relatively moist 
(humid) air which is entrained and lofted from near the 
Earth's surface to higher and colder regions of the 
atmosphere. 

Particle is a generic term which includes dust, pebbles, rocks, 
boulders, water, ice, or any other condensed phase mass in the 
cloud. The lofted particles are of many different sizes and shapes. 
For convenience, the size of a particle is usually taken as the 


equivalent spherical diameter, a: 


a 6m \N3 Diameter (m = 1 gm when a = 1 cm and 
Oo e. = 1.9 gm/cm”) 
where 
™ Mass of particle 
OG Density of particle 


*Mass concentration is also commonly called density. 


The weapon yield (W) is the energy released in a nuclear 


detonation.* The yield is usually expressed in equivaient tons 


of TNT, where 


lL ton of TNT is assumed to 
4.2 x 10° Joules 


1 kiloton 


release 10° calories or 


L kt = 4.2 x 10!?7 Joules 


1 megaton = 1 Mt = 4.2 x 10!5 Joules 


The first nuclear detonation, shot TRINITY in 1945, had a yield 
of ~19 kt. The yields of the weapons detonated over Hiroshima 


and Nagasaki were approximately 13 kt and 23 kt, respectively. 
Shot CASTLE BRAVO in 1954 at Bikini (in the Pacific Proving Ground) 


had the largest U.S. test yield 


of ~15 Mt. 


The height of burst (HOB) is the device height above the 
Earth's surface. The scaled height of burst (SHOB) is defined 


as 


SHOB = 


Buried bursts have negative HOB 


A contact surface burst is 
A nuclear detonation this close 
dynamic cratering process which 


material into the atmosphere. 


HoB /w!/3 (2-1) 


and SHOB. 


defined as |SHOB] < 5 ft/kt'/?. 
to the ground surface causes a 


ejects large amounts of earth 


A contact surface burst couples 


about 8% of the nuclear yield to the ground in a few microseconds. 


The details of the energy coupled to the ground, and the amount 


of ejecta thrown into the atmosphere will depend on W, HOB, the 


device construction, and the properties of the ground material. 


*The nature of the nuclear energy released and a general descrip- 
tion of nuclear weapons effects are in Reference 1. 


As the SHOB increases above 5 ft/kt (1.6@., about 50 ft. for 
W = 1 megaton), the fraction of energy coupled directly to the 
ground and the amount of ejected mass begins to decrease rapidly. 
An estimate of the crater ejecta mass source vs HOB and SHOB for 
W = 1 megaton is shown in Figure 1 below. All other sources of 
ground material which become entrained in the nuclear cloud come 
from near the ground surface and can be conveniently combined into 
the sweep-up mass source, The sweep-up mass source dominates 
the lofted cloud mass characteristics for SHOB 2 20 ft/kt™’. 
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Fig. 1. Crater Ejecta Mass Source with Ejection Speeds of Over 
10 m/s versus HOB for W = 1 Mt.? 


SECTION 2 


DUST/DEBRIS SOURCES 


2.1 CRATER EJECTA MASS SOURCE 


For contact surface bursts, crater ejecta represents the 
primary source of ground material which becomes entrained in the 
nuclear cloud. The following figure indicates the general physical 
phenomena, pressures, and times involved in crater formation for a 


nominal 1 megaton contact surface detonation. 
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Fig. 2. Typical Pressure-Time Sequence for Energy Coupling and 
Cratering. 
The following photos show the fireball and ejecta from the 
Slightiy buried burst JOHNIE BOY (W = 0.5 kt, HOB ~-2.0 ft, at 
NTS). Note that some of the crater ejecta emerges from the 


early-time fireball in distinct rays. 


Interacting 
with Ground 


0.2 sec. 


igs 3. Photographic Exampte for 


Shot, JOHNIE BOY (W = 0.5 kt, PLATE 14 


HORIZONTAL SCALE 5,7 meters/mm 
DOB = .6 m).” 


iz 


0.8 sec. 


2.0 sec. 4,0 sec. 


SCALE (.8-4 sec.) 4.5 meters/mm 


8.0 sec. ; 
4. Photographic Example for 12. sec. 20. sec. 
Shot JOHNIE BOY (W = 0.5 kt, PLATE 15 SCALE ( 8-20 sec.) 8.5 meters/mm 
DOB = .6 m).* 


Theoretical calculations (or computer code simulations) of 
the crater and ejecta formation use finite difference analogues 
of the differential equations for conservation of mass, con- 
servation of momenta, conservation of energy, and approximate 
equations of state for the weapon and ground materials. 


These numerical simulations predict the particle velocity, 
gensity and temperature of the incipient ejecta versus radius and 
time. (Most crater/ejecta calculations are axisymmetric and ejecta 
rays are ignored.) Figure 5 (next page) shows a velocity vector 
field near the dynamic crater at t ~0.2 seconds for a 1 Mt contact 
burst over the indicated geology. The incipient ejecta is flowing 
through the original ground surface (Z = 0) at radii between 40 m 
and 80 m at 0.2 seconds. The ejection angles are roughly 45°. 


The ejecta mass rate per unit area flowing across Z = O at any 
radius is pV. At t = 0.2 seconds, pV peaks at a radius of about 
50 meters and has an associated vertical velocity of ~60 m/s. 
Figure 6 shows the 


analogous vertical 


velocity (i-.e., Cratering 
Calculations 


associated with 


the peak mass 


flux) as a func- ,"° 

tion of time : 

from 5 nuclear = Velocity x 2 
cratering cal- re Vertical Yelocity(¥,,] 


culations. The 
various calcul- 
ations show 
variations of 


about a factor 


of 2 from the oe “mt ; 
TIME (sec) 
mean curve. 
Fig. 6. Ejecta Verticdl Velocity at a Range of 
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Fig. 5. Crater Ejecta Source Example: Velocity Field at t ~0.2 Seconds for 
a W = 1 Mt Surface Burst. ® 


Of particular importance (for the lofted mass characteris-— 
tics of the nuclear cloud) is the particle size distribution of 
the incipient ejecta mass. This size distribution is not provided 
by current cratering calculations. Physically, the particle size 
distribution will depend on, at least, the stress, strain, and 
temperature history of each element of the incipient ejecta. In 
general, as the ejecta flow into the air, the particle size dis- 
tribution will vary with both the time of ejection and the range 
from the burst point. The incipient ejecta initially consist of 
vaporous soil and bomb debris material. At somewhat later times, 
melted soil flows into the fireball. The bulk of the ejecta, 
however, is composed of solid material which has been shocked to 
some peak pressure, distorted, and flowed by the cratering dynamics. 
For example, about 500 kton of soil experiences 7 kbar (100 ksi) or 
more of pressure in a 1 Mt yield surface burst. Thus, the crater 
ejecta particle size uncertainty involves not only natural in-situ 
size variations, but also uncertainties in soil/rock break-up and 


agglomeration during crater formation. 


The crater does not completely form until about 5 seconds for 
a 1 Mt surface burst (an example of a predicted final crater pro- 
file is shown with a dashed line on Figure 5). Roughly 2 Mtons of 
ejecta mass (1 Mton = 10‘? gm) is thrown from the crater, however, 
note that only a small fraction* of the ejecta mass becomes lofted 
as part of the nuclear cloud. Most of the mass falls back to earth 


between 10 seconds and 1 minute. 


*Jt is estimated that ~0.3 Mtons of ejecta mass per Mt of weapon 
yield remains aloft at 5 minutes after detonation of a contact 
surface burst. 


2.2 SWEEP-UP MASS SOURCE 


For non~cratering HOB* detonations, the dust/pebbles which 
become entrained in the nuclear cloud originate very near the 
ground surface (probably within a few centimeters or so). A 
"dusty" sweep-up layer is formed as this near-surface material 
mixes with the high velocity winds behind the air shock wave. 
The sweep-up layer is fairly persistent, and some of the dust/pebble 
particles in this layer are subsequently drawn up into the main 


cloud. 


Shock Wave 


Thermal Pre-Shock Layer 
+? 


Thermal Precursor 


{ 


Thermal heating of an air layer near the ground surface may 
occur prior to shock arrival; this heated layer causes a thermal 
precursor flow field which appears to enhance the sweep-up layer. 
A thermal precursor is a pressure and velocity wave disturbance 
which travels faster in the heated near-surface thermal layer 


than in the cooler air above the surface. 


*HOB is commonly used as an abbreviation for "height of burst" as 
well as the burst distance above the ground surface. 


Figure 7 is a photographic example of a relatively late-time 
(~4 sec) sweep-up layer for shot GRABLE (W = 15 kt, HOB = 524 ft, 
SHOB = 212 ft/kt'#). Based on existing Nevada Test Site (NTS) data 
for SHOB between 150-500 ft/kt!, the maximum radial extent of the 
Ssweep-up layer is about aa ~1000-1300 Ft/kt 8 which corresponds 
to a peak overpressure of ~5-10 psi. The maximum height of the 
sweep-up layer varies with the SHOB and ground surface character- 
istics, but a few hundred feet are typical for the relatively 
tow yield NTS shots. Techniques for scaling the maximum radius and 
height of sweep-up layers to larger yields (say W 2 1 Mt) are still 
uncertain; however, estimates of ~1000 ft heights for W = 1 Mt seem 
plausible. 


The early-time growth of the sweep-up layer is clearly coupled 
to the shock wave and air flow characteristics near the ground 
surface. Figure 8, drawn from photographs at 3 times after burst, 
illustrates the fact that the sweep~-up layer begins forming very 
quickly behind the shock wave. 
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Fig. 8. Tracings from Technical Photography: Shot GRABLE.’ 
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The mass concentration or density (9) of the dust which 
becomes entrained in the sweep-up layer has been measured in two 
nuclear tests. For example, a density of about 2 x 107? gm/cm’ 
has been measured in shot TEAPOT MET (W = 22 kt, SHOB = 140 f£t/kt1/*) 
at a range corresponding to peak overpressures and particle 
velocities of ~20 psi and ~600 ft/sec, respectively. This density 
was measured at a 3 ft height above the ground surface. Uncer- 
tainty factors of about 3 are suspected in this dust density 


measurement. 


The actual concentration profile above the ground surface is 
not known: however, the following curve shows the general order 
of magnitude variation expected near the ground surface at a 


typical range and time in the sweep-up layer. 
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The high dust density region near the original ground surface 
represents the mass source for the sweep-up layer and ultimately 
the dust source for the nuclear cloud. The general nature of the 
interactions between the air and the near-surface material are 
sketched below: 


Air Velocity 
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The multiphase physics (air plus solid/liquid soil particles) 
in the first 1-10 cm above and below the original ground surface 
involve viscous stresses which are associated with high gradients 
in horizontal air velocity near the ground surface. Also, the soil 
permeability and strength properties will determine the transfer 


of alr momentum and energy to the near-surface soil material. 


Immediately behind the shock wave on the ground, there will 
be some permeation (i.e., diffusion) of air into the porous sur- 
face material. The near-surface (~l cm) permeability variations 
with depth and range will lead to scouring and lofting of near- 
surface material on a short-time scale. On a longer time scale, 
air will permeate through the pores deeper into the soil. As the 
overpressure on the surface subsequently decays, the air pressure 
in the pores below the surface will be higher for a time than the 
pressure being applied onto the surface by the decaying airblast 
wave. When the air pressures within the pores of the soil exceed 
the applied overpressure, there is an unbalanced upward force on 
material near the surface. ‘Thus, soil material is ejected upward 
into the near-surface air flow and sweep-up layer. This effect 
has been observed in laboratory experiments and in HE field 


experiments, and has been theoretically calculated. 


The deeper pore-air flow phenomenon may eject ~10 cm of soil 
material into the sweep-up layer.® There is probably a sufficient 
delay behind the shock wave that this material will not be 
accelerated to high velocities for a single burst. On the other 
hand, in a multiple nuclear burst environment, this phenomenon 
May represent the most significant dust/pebble/rock source for 
subsequent bursts at some sites. The permeability and cohesive 
strength of the shocked and heated ground material will probably 
be the most important site dependent physical characteristic for 


this phenomenon. 


Once the dust/pebbles become lofted above a few meters, the 
complex air flow field behind the shock wave will dominate the 
further rise and turbulent diffusion of the dust/pebble mass. 


The air flow field behind the near-surface shock wave is 
quite complicated due to 


@ unsteady pressure and flow dynamics including Mach stem 
and thermal precursor formation 


@e turbulent boundary layer effects involving in-situ 
bushes, trees, rocks/boulders, and smali~-scale 
hills/valleys 
Many questions and issues remain unanswered in predicting the 
formation and nature of the sweep-up layer; Figure 10 summarizes 


some of these issues. 
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Fig. 10. Issues Concerning Sweep-Up Layer Dust and Pebbles. 


SECTION 3 


NUCLEAR CLOUD FORMATION 


3.1 EARLY-TIME OR FIREBALL PHASE (t * 10 seconds for W = 1 Mt) 


The very early-time phase (t * 0.1 second) of a nuclear 
cloud is characterized by extremely high fireball temperatures 
(10°-10°°K) and shock wave pressures (107-10° psi) as indicated 
in the following Figures lla and llb. These curves are from a 
numerical simulation of a 2 Mt free air burst; this burst is 
equivalent to a 1 Mt contact surface burst over an idealized rigid 


surface with no crater and no ejecta. 
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The late stages of the fireball growth and associated shock 
front behavior are indicated on Figures llc to lle. The per- 
sistence of the high fireball temperatures and low air densities 
(Figure 1ld) set the stage for the buoyancy dominated rise of the 
hot atr and the lofting of entrained dust/debris. 
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Early-Time Contact Surface Bursts 


The detonation of a nuclear device very near the Earth's 
surface causes ground material to be thrown from the developing 
crater into the fireball. The initial pressures in the ground 
start at about 1000 Mbars (~1.5 x 10'° psi),* and consequently a 
strong shock wave develops and propagates into the local geologic 
material as described in Section 2.1. The geology is, therefore, 


quite important in determining the crater ejecta characteristics. 


The ejecta flows into the fireball forming a dynamic multi- 
phase region involving aerodynamic drag and thermal interactions. 
The thermal interactions include vaporization of solid/liquid 
material. The vaporization (and later recondensation) of the 
earth material is potentially significant when comparing shots 


on land with shots on water. 


Figures 12a and 12b show predicted air velocities and temp- 
eratures 1 second after a 1 Mt contact surface burst on a soil/rock 
geology; estimated models were used for crater ejecta character- 
istics and multiphase drag/thermal interactions. Ejecta mass 
groups representing different size particles were assumed in the 
numerical calculations and each of these size groups will inter- 
act differently with the air. The boundaries of various size 
groups are shown on Figure 12a with the larger size particles 
traveling the larger distances from the burst point. The ejecta 
influences both the air velocities and temperatures. The temp- 
eratures (Figure 12b) are reduced in the interior of the fireball 
by relatively cold ejecta which mixes with the hotter air. 
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dirt concentration field 
at t = 14 seconds. The 
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Early-Time HOB Detonations 


As the height of the nuclear burst is increased above the 
Earth's surface, less of the weapon's energy is coupled to ground 
material and crater ejecta begins to play a smaller and smaller 
role in forming the nuclear cloud. This transitional region will 
not be discussed in detail; its characteristics are a combination 
of the contact surface burst phenomena discussed previously 
and the non-cratering HOB physical phenomena described in the 


remainder of this subsection. 


When the spherical air shock wave from the burst interacts 
with the Earth's surface, complex reflectéd waves are transmitted 
back into the previously shocked air. These reflected waves 
strongly influence the subsequent flow field of the air and lofted 
material. Figures 13a and 13b show an example for shot U/K GRABLE 
(W = 15 kt, SHOB ~200 ft/kt!*) of predicted air velocity vectors 
just prior to and shortly after shock reflection off the ground. 
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3.2 CLOUD RISE AND STABILIZATION PHASES (10 sec S$ t $ 10 min) 


After approximately 10-15 seconds for a 1 Mt burst, buoyancy 
forces begin to dominate the air velocity field characteristics. 
Buoyancy forces become important because the atmospheric pressure 
field returns to nearly ambient conditions, but peak air tempera~ 
tures still exceed 1000°K at 10 seconds. The low density fireball 
air will accelerate upward causing a vortex air flow field to 
develop which entrains relatively cool air, earth material, radio- 
active material, and atmospheric humidity. The vortex flow carries 
the entrained material through the stem and into the upper portions 
of the cloud. 
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At any time after burst, the amount of earth material and 
the spatial distribution of this material in the cloud depend on 
the weapon yield, HOB, and material properties (e.g., size dis- 
tribution and strength) of the soil/rock. An example follows: 


Figure 15 is an estimated spatial distribution of the lofted 
dust and pebble mass for a 1 Mt contact surface burst at 5 minutes 
after detonation. Figure 15a shows the mass concentration field. 
The altitude of the cloud is roughly 70,000 ft with a radius of 
approximately 20,000 ft at 5 minutes. Note that a vortex region 
of relatively high density (107° gm/cm*) develops at an altitude 
of about 50,000 ft and a radius of between 10-15,000 £t. 


The nature of the particle size distribution of this mass is 
presented in Figure 15b, 


which shows the cumulative & bust (D,6.5 m) 


O Pebbles (,5 em¥ D,*5 em) 
@ Dust and Pebbles 


lofted mass above any alti- 
tide Z at 5 minutes. Note 
that most of the mass is in 
the main cloud {about 210 
kilotons of lofted mass above 
40,000 ft}; the stem contains 
about 50 kilotons of mass. 
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1a. lofted Mass Conceatration 15b. Cumulative Dust and Pebble Mass 
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larger particles have Fig. 15. Lofted Dust/Pebble Spatial (R,Z) 
Distribution for a 1 Mt Surface Burst 
landed on the ground by at 5 minutes. > 


5 minutes: 


The buoyancy dominated phase ends when the cloud reaches its 
maximum or stabilized cloud altitude. This time is typically 3 to 
53 minutes. While the cloud height will stop growing, and in fact 
decrease in many cases after a few minutes, the cloud diameter 
will continue to increase; 4 20 m/sec 
the cloud diameter will 
increase until about 7-10 
minutes even in the absence 


of atmospheric winds. The 


Vortex 


increase in cloud diameter 
is due to the post- 
stabilization air flow 
characteristics which occur 
near the cloud top as 


illustrated in Figure lo. 


The reverse vortex forms <A, 
Primary 


as a result of downflow 
of relatively dense air 


above the cloud top. 
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Relatively dense air from 
lower altitudes is "pushed" 
upwards to lower density 
portions of the atmosphere; 
when the upward velocities 
caused by buoyancy get 
sufficiently small, gravi- 


tational forces on the 


relatively dense air 


dominate the flow field. 40 
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Fig. 16. Air Flow Field at 5 minutes 
for a 1 Mt Contact Surface Burst.’° 


3.3 EXPERIMENTAL CLOUD DIMENSION DATA 


Figure 17 is an example of nuclear cloud dimensions versus 
time for shot CASTLE BRAVO (W = 15 Mt surface burst in the Pacific). 
At about 5 minutes, the cloud top reaches a maximum altitude of 
114 kft; at 9 minutes, the cloud and stem diameter are about 


330 kft and 28 kft, respectively. 


Figure 18 shows experimental data for cloud altitude versus 
weapon yield at 1 minute after detonation; Figure 19 shows the 
maximum (or stabilized) cloud altitude versus yield. The various 
symbols on these figures indicate whether the shot was conducted 
in the Pacific or at the Nevada Test Site, and the type of 
measurement technique used. The scatter in the data is due pri- 
Marily te variations in atmospheric conditions (e.g., winds, 


lapse rate, humidity) and/or measurement errors. 


Figures 20 and 21 show experimental data for cloud diameter 
versus weapon yield at 1 minute and 10 minutes after detonation, 
respectively. Atmospheric winds and turbulence will often cause 


continued growth in the cloud diameter after 10 minutes. 
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Fig. 18. Cloud Height Above Burst Point Versus Weapon Yield at 
t = 1 Minute. ?? 
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Fig. 19. Cloud Maximum/Stabilization Height Above Burst Point 
Versus Weapon Yield. ** 
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Cloud Diameter Versus Weapon Yield at t = 1 Minute.?? 
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Fig. 21. Cloud Diameter Versus Weapon Yield at t = 10 Minutes.’? 
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